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In order to study the pairing symmetry in non-centrosymmetric superconductors, we solve the 
linearized Eliashberg's equation on the two-dimensional extended Hubbard model in the presence 
of the Rashba-type spin-orbit coupling (RSOC) within the random phase approximation. In the 
presence of the RSOC, three types of pairing symmetries appear in the phase diagram with respect 
to the on-site Coulomb repulsion U and off-site one V. Each of pairing symmetries is admixture of 
spin-singlet and -triplet ones. On the basis of analytical study, it is found that the admixture of 
spin-singlet and -triplet components depends on not only the predominant pairing symmetry but 
also dispersion relation and pairing interaction. 



I. INTRODUCTION 

Since the discovery of superconductors without the 
inversion symmetry in CePt 3 SiJ2 non-centrosymmetric 
superconductors has been studied intensively as un- 
conventional superconductors. In particular, theo- 
retical studies suggest interesting properties of non- 
centrosymmetric superconductors such as the magneto- 
electric effect pHSl anisotropic spin susceptibi lity,^ ^ anc j 
the helical vortex state in magnetic fields! 12 * 13 * Today, 
there are various relevant systems in two-dimensional 
electron gas at heterointerface, e.g. between SrTiOa 
and LaA10 3 P and non-centrosymmetric crystals, e.g. 
CePt 3 SiP UIrP CeRhSia, 111 CelrSiaF 1 CeCoGeaf™ 
and L^Pd.Pta-.BpHU] 

In superconductors with the inversion symmetry, pair- 
ing symmetry is classified into even- and odd-parity, i.e. 
spin-singlet and -triplet. In the non-centrosymmetric su- 
perconductors, on the other hand, a dmixture of spin- 
singlet and -triplet pairings is realized? 2 * 7 * 8 * 23 * 24 ! The ad- 
mixture is induced by the antisymmetric spin-orbit cou- 
pling, which is generated by the lack of the inversion sym- 
metry. For example, it has been proposed that admix- 
ture of spin-singlet s-wave and spin-triplet p-wave pair- 
ings is realized in a non-centrosymmetric he avy ferm ion 
superconduc tor C ePtaSi from both theoretica l 11 * 25 * 26 -* and 
experimental 1 "* 2 ^^ studies. 

In particular, non-centrosymmetric heavy fermion 
superconductors, e.g. CePtaSiP UIrP CeRhSia! 111 
CelrSiaJ^ and CeCoGe 3 , 1 18 * 19 * are of interest because 
the superconductivity originates from the Coulomb re- 
pulsion. However, there are few theoretical studies 
on these materials on the basis of the microscopic 
calculation! 11 * 25 * 3 - 3 -* 3 ^* It is desired to microscopically 
understand the admixture of spin-singlet and -triplet 
pairings in non-centrosymmetric superconductors where 
Coulomb repulsion mediates pairing in more detail. 

In order to study the above admixture, it is valuable to 
employ the extended Hubbard model because this model 
exhibits various pairin g sym metries in the presence of 
the inversion symmetry! 35 " 47 * On the extended Hubbard 
model, the off-site Coulomb repulsion is considered in ad- 
dition to the on-site one. It is well known that, while the 



on-site Coulomb repulsion induces the spin fluctuation, 
the charge fluctuation coexists with the spin fluctuation 
by introducing the off-site Coulomb repulsionP^SH j} ue 
to the coexistence of the charge fluctuation with the spin 
one, especially on a two-dimensional square lattice near 
half-filling, three types of pairing symmetries, i.e. spin- 
singlet c? a; 2_j ; 2-wave, spin-triplet /-wave, and spin-singlet 



ci^-wave ones, compete against each others! 3 - 5 -* 

In the present study, in order to clarify pairing symme- 
try in non-centrosymmetric superconductors where the 
Coulomb repulsion mediates pairing, we investigate the 
two-dimensional extended Hubbard model in the pres- 
ence of the Rashba-type spin-orbit coupling (RSOCJpHl 
on the basis of the random phase approximation (RPA) . 
The RSOC induces breakdown of the inversion symmetry 
and admixture of pairing symmetry. 

The admixture of pairing symmetry has already been 
studied on the extended Hubbard model in the presence 
of the RSOC on the basis of the RPA by Yokoyama et 
ali^ However, they investigated only the region where 
the off-site Coulomb repulsion is small, i.e. spin-singlet 
d x 2_ y 2-w&ve pairing state is stable in the absence of the 
RSOC. Moreover, there were two simplifications in the 
pairing interaction. One is that they neglect cross terms 
of the bubble- and ladder-type diagrams, which are gen- 
erated by the off-site Coulomb repulsion and the RSOC. 
The other is that the ladder-type diagrams with the off- 
site Coulomb repulsion are excluded. In the present 
study, there is no simplification described above. 

The present paper is organized as follows. In SjTTJ we 
formulate the linearized Eliashberg's equation on the ex- 
tended Hubbard model in the presence of the RSOC on 



the basis of the RPA. In i III A we show results obtained 



by the numerical calculation within the RPA. After that, 
we discuss the pairing symmetry on the basis of analyti- 
cal study in §HIB| The summary is given in §IV| 



II. FORMULATION 

We start with the two-dimensional extended Hubbard 
model in the presence of the RSOC. The Hamiltonian is 
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given by 

W =^£kci s c ks + U^n^riii + V ^ n ia nj a > 

k.s i (i,j),s,s' 

-Xj2l9(k)-B-] ss/ cl s c ks ,, (1) 

k,s.s' 

Ek = — 2t(cos k x + cos k y ) — At' cos k x cos k y — fi, (2) 

where is an annihilation (a creation) operator for 
an electron with spin s and momentum k, rii s is a 
number operator for an electron with spin s at site i, 
and denotes a set of the nearest neighbor sites. 

Ek is the dispersion relation, where v- 1 ' is the (sec- 
ond) nearest neighbor hopping on a square lattice and 
fi is the chemical potential. We consider the on-site 
Coulomb repulsion U and the off-site one V between 
the nearest neighbor sites. The fourth term is the 
RSOC, where A is a magnitude of the RSOC and <x 
are the Pauli matrices. The vector g(k) with the re- 
lation g(k) = —g(—k) induces breakdown of the inver- 
sion symmetry. We adopt g{k) = (~v y (k),v x (k),0)/v 
with the quasiparticlc velocity v x ( y )(k) = d£k/dk x ( y ) — 

g{k) is normalized 



2t sin k„ 



(y) 



At' Sin k X (y) COS ky( X ). 



by the average velocity v which is given by v 2 — 
J2k \_ v x{k) 2 + Vy(k) 2 ] /N, where TV is the number of k- 
meshes. The bare Green's function is given by the fol- 
lowing 2x2 matrix in spin space, 



G(k) 



G n (k) 



G U (k) 



(iw„ - £ fe )I + \g(k) ■ <t 



(3) 



where I is a unit matrix and k = (iw„, fe) is an abbrevi- 
ation. u> n — (2n — l)irT is the Matsubara frequency for 
fermions, where n is an integer and T is temperature. 

In order to estimate the pairing instability, we solve 
the linearized Eliashberg's equation within the RPA 

aA slS2 (k) = -^ E r s;? a .,.<(*-*') 

fc / ,Cl,C2, S%,S4 

xP^(k')P C2 (-k)F S3Si (k'), (A) 
F S1S2 (k) = £ G S1S3 (k)G S2Si (-k) A S3S4 (fc), (5) 

S3,S4 



restarted Arnoldi method.^ This method is powerful in 
solving an eigenvalue equation with nearly degenerate so- 
lutions. 

The spin-singlet and -triplet components with S z = 
are extracted by [A^(fc) ± A^(fc)] /2, where the sign 
+(— ) corresponds to spin-triplet (-singlet) one. The spin- 
triplet components with S z — ±1 are given by A^m.) (&) 
for the sign +( — ). In the present paper, we choose the 
solution whose spin-singlet component is real. The spin- 
triplet (S z — ±1) components are imaginary. The real 
and imaginary parts of the spin-triplet (S z = ±1) com- 
ponents have same amplitude while nodes of the real 
and imaginary parts have the relation of rotation around 
fe=(0,0). 

Within the RPA, the effective pairing interaction 
PsjS2s 3 s 4 (9)' wnere 1 = (i^mj?) with the Matsubara fre- 
quency for bosons v m = 2mirT, is obtained by collecting 
the infinite series which consist of the irreducible suscep- 
tibility in the diagrammatic expression. The irreducible 
susceptibility is given by 



^o!si 2 s 2 s 3 s 4 ( ( J l ) 



^E Gs ^(fc + <7)G S4S2 (fc) 

k 



x P Cl (k)P C2 (k). (6) 



The dressed susceptibility is given by 



x{q)=xo(q) i - f o(<z)xo(<?) 



(7) 



where the matrices are 20 x 20 ones with spin indices Sj 
and phase factor ones Cj defined as 



M 



M c 



( M 11 

M 21 
M 31 
M 41 

Vm 51 

M ClC2 
\ M ClC2 



M 12 
M 22 
M 32 
M 42 
M 52 



M 13 
M 23 
M 33 
M 43 
M 53 



M 14 
M 24 
M 34 
M 44 
M 54 



M ClC2 
M ClC2 
M ClC2 
M Cl ° 2 



M': 



C1C2 

M ClC2 
M tUt 

M ClC2 
M Cl ° 2 



M 15 \ 

M 25 
M 35 
M 45 
M 55 J 

M ClC2 
iw ttW 

M ClC2 
M ClC2 
M ClC2 

■\r-\r-\r-\r 



(8) 



(9) 



The matrix r (q) is given by 



where P(k) = (1, cos k x , sin k x , cos k y , sin k y ) is the phase 
factor which originates from ladder-type connections 
with the off-site Coulomb repulsion V(q) — 2X^(008^ + 
cosq y ) in the diagrammatic expression. The linearized 
Eliashberg's equation Q and ^ is an eigenvalue equa- 
tion whose eigenvalue and eigenfunction are a and A(fc), 
respectively. When the eigenvalue a reaches unity, tem- 
perature T corresponds to the superconducting transition 
temperature To Thus, the eigenvalue a implies the pair- 
ing instability with the gap function A(fc). In solving the 
linearized Eliashberg's equation, we employ the implicit 



T^ 2 (q) 



n iC2 (q) 




(ci=C2=2-5), 
(ci # c 2 ). 



(10) 



(11) 
(12) 



3 



By using the dressed susceptibility, the effective pairing 
interaction is expressed as 



A S1S2S3S 



T (q)x(q)T {q) 
where the matrix T' (q) is given by 



cic 2 

S1S3S4S2 



ylc 1 c 2 



(13) 



-V(q) 

-U-V(q) 

-U-V(q) 






f' ClC3 (q)=6 (d^l.ca^l). 

Fig. [T] shows a sense of the formulation for the pairing 
interaction within the RPA in a diagrammatic represen- 
tation. In the present RPA for the extended Hubbard 

FIG. 1: Diagrammatic sense of the effective pairing interac- 
tion r(g) within the RPA. To(q) includes the on- and off-site 
Coulomb repulsions. 



model, cross terms of the bubble- and ladder-type dia- 
grams, which generated by the off-site Coulomb repulsion 
and the RSOC, are taken into account. 

Within the RPA, the spin susceptibility is expressed as 



X%(q) =^J dre i ^(T T S^r,q)S'\-q)) 

a / j u s 3 s 4 As 1 s 2 s 3 s i \H/' J s 2 Si' 



where 



S*(T,q) =e Hr St(q)e- nr , 



(16) 
(17) 

(18) 
(19) 



with £,?y = x,y,z. Similarly, the charge susceptibility is 
expressed as 



Xch(q) 



1 

~2N 



dTe i ^(T T p(r,q)p(-q)) 



' n ^ ] XsSS's' (?) ' 



where 



p(r, q) =e HT p(q)e 



-Ht 



(20) 
(21) 

(22) 
(23) 



In the present paper, we choose t — 1 for a unit of 
energy. The second nearest neighbor hopping, tempera- 
ture, and filling are always t' — 0.1, T — 0.04, and 0.8 
electrons per site, respectively, in the actual numerical 
calculation. We take 64 x 64 fc-meshes and 1024 Matsub- 
ara frequencies. 



III. RESULTS 

A. Numerical calculation within the RPA 

In this subsection, we show results obtained by the 
numerical calculation within the RPA. 



1. In the absence of the RSOC 

First, we check the pairing symmetry in the absence of 
the RSOC, i.e. A = 0. Fig. [2] shows U-V phase diagram. 
We identify the pairing symmetry with the largest eigen- 




FIG. 2: (Color online) U-V phase diagram in the absence 
of the RSOC within the RPA. The broken line shows V = 
-([/ - 2.8)/4 between points A (2.8,0) and B (0.9,0.475). 



value a in the linearized Eliashberg's equation (|4|) and 
([5]) at T — 0.04. The boundary with the spin- or charge- 
density-wave (SDW or CDW) phase is determined by the 
Stoner's factor, which is defined as the largest eigenvalue 
of the matrix f o(q)xo(q) in the dressed susceptibility Q. 
When the Stoner's factor reaches unity, the dressed sus- 
ceptibility diverges. In the present paper, we define the 
SDW or CDW phase as the region where the Stoner's fac- 
tor reaches 0.98. As shown in Fig. [2j three types of the 
pairing symmetries can appear by tuning U and V; spin- 
singlet c? a; 2_ l) 2-wave, spin-triplet /-wave, and spin-singlet 
d xy -wsve pairing symmetries. 

Three types of the pairing symmetries are caused by 
the cooperative/competitive spin and charge susceptibil- 
ities controlled by U and Fig. [3] shows the spin 
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and charge susceptibilities in the absence of the RSOC 
on the broken line V = -{U - 2.8)/4 in Fig. § (U-V 
space). Here, we plot the maximum values of Xgpvl) an d 




B U A 



FIG. 3: (Color online) The spin and charge susceptibilities 
on the broken line V = -(U - 2.8)/4 in Fig. [5] (U-V space) 
without the RSOC. The maximum values are plotted for each 
(U, V)-point. The result is numerically obtained within the 
RPA. Points A and B correspond to those in Fig. [2] 

Xch(q) for each (U, V^-point. Note that, in the absence 
of the RSOC, = x»(g) = X z s ;(q\ and yto(q) = 

for t; =/= T], i.e. the spin susceptibility is isotropic. In the 
region where U (V) is large (small) on the broken line 
in Fig. [2j the spin susceptibility is dominant as com- 
pared to the charge one. In the region where U and V 
are intermediate on the broken line in Fig. [2j the spin 
and charge susceptibilities are comparable. In the region 
where U (V) is small (large) on the broken line in Fig. 
[2j the charge susceptibility is dominant as compared to 
the spin one. Thus, surveying from point A to B along 
the broken line in Fig. [2] dominant fluctuation changes 
from spin one to charge one. 

In the region where U (V) is large (small), dominant 
spin susceptibility, which is repulsive for spin-singlet pair- 
ing channel, generates spin-singlet d x 2_ y 2-wave pairing 
symmetry. This pairing symmetry has the sign change on 
the Fermi surface through the nesting vector, where the 
spin susceptibility has peak in momentum space. Thus, 
for spin-singlet d x 2_ y 2-vr&ve pairing symmetry, the spin 
susceptibility works attractively. In the region where U 
and V are intermediate, spin-triplet pairing channel ex- 
ceeds spin-singlet one in the pairing interaction and spin- 
triplet /-wave pairing symmetry is favored. The reason 
why spin-triplet pairing channel exceeds spin-singlet one 
is as follows. In spin-singlet pairing channel, the spin 
and charge susceptibilities are competitive in the pairing 
interaction. In spin-triplet pairing channel, on the other 
hand, the spin and charge susceptibilities are cooperative. 
In addition, since the pairing interaction for spin-triplet 
channel is originally attractive, stable pairing symmetry 
is /-wave one, which has no sign change on the Fermi 



surface through the nesting vector, where the spin and 
charge susceptibilities have peaks in momentum space. 
In the region where U (V) is small (large), spin-singlet 
pairing channel and spin-triplet one are comparable in 
the pairing interaction. Then, the charge susceptibility is 
attractive for both pairing channels and spin-singlet d xy - 
wave pairing symmetry is favored. This pairing symme- 
try has no sign change on the Fermi surface through the 
nesting vector, where the charge susceptibility has peak 
in momentum space. The total number of nodes in the 
gap function for spin-singlet d xy -wave pairing symmetry 
is less than that for spin-triplet /-wave one. Details have 
been discussed by Onari et aiP2 



2. In the presence of the RSOC 

Next, we introduce the RSOC. It has been known that 
the RSOC makes the spin susceptibility anisotropic tj-e. 
X^(q) ^ X v 4{q) ¥= Xi;{q) and X f£(<?) ^ for ^ ^ In 
the presence of the RSOC (A = 0.3), the pairing sym- 
metry is shown in Fig. [4j The pairing symmetry and 




FIG. 4: (Color online) U-V phase diagram in the presence of 
the RSOC (A = 0.3) within the RPA. The broken line shows 
V = -(C/-2.8)/4 between points A (2.8, 0) and B (0.9, 0.475). 

boundary with the SDW or CDW phase are determined 
in the same manner as in the absence of the RSOC. As 
shown in Fig. [4j three types of pairing symmetries can 
appear by tuning U and V. In the region where U (V) 
is large (small), predominantly spin-singlet dj,2_ y 2-wave 
pairing symmetry admixed with spin-triplet (S z = ±1) 
p-wave one is the most stable. We call this pairing sym- 
metry p + D x 2_ y 2--wsve one in the present paper. Mo- 
mentum dependence of the gap function is shown in Fig. 
[5] (a). In the region where U and V are intermediate, pre- 
dominantly spin-triplet (S z = ±1) /-wave pairing sym- 
metry admixed with spin-singlet c? a; 2_ l) 2-wave one is the 
most stable. We call this pairing symmetry d x i_ y i + F- 
wave one in the present paper. Momentum dependence 
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FIG. 5: (Color online) Momentum dependence of the gap 
functions at u) n — irT in the presence of the RSOC (A = 0.3) 
within the RPA. (a) p + D x 2 _ y n -wave pairing symmetry for 
(U,V) = (2.8,0), (b) d x 2_ y 2 + F-wave one for (U, V) = 
(2.1, 0.6), and (c) p+D^-wave one for (U, V) = (0.8, 0.45) are 
exhibited. Panels on the left (right) side show the predom- 
inant components (subcomponents). Black solid lines and 
green broken ones denote the Fermi surface and node of the 
gap functions, respectively. For spin-triplet components, only 
ReA-|-|-(fe) is shown. 



of the gap function is shown in Fig. [5] (b). In the re- 
gion where U (V) is small (large), predominantly spin- 
singlet d xy -wave pairing symmetry admixed with spin- 
triplet (S z — ±1) p-wave one is the most stable. We call 
this pairing symmetry p+D xy -wave one in the present pa- 
per. Momentum dependence of the gap function is shown 
in Fig. [5](c). Focusing on the predominant components, 



the pairing symmetries are the same as in the absence of 
the RSOC, i.e. spin-singlet d x 2_ y 2-wave, spin-triplet /- 
wave, and spin-singlet d xy -w&ve pairing symmetries, ex- 
cept lack of spin-triplet component with S z = 0. Phase 
boundaries are almost unchanged by the RSOC as shown 
in Figs. [2] and [2] Namely, in the presence of the RSOC, 
subcomponents are admixed with predominant compo- 
nents whose pairing symmetries are determined by coop- 
erative/competitive spin and charge fluctuations without 
the RSOC as discussed above. It is discussed later in 
pil B| how the pairing symmetries of the subcomponents 
are determined. 

The eigenvalues a in the linearized Eliashberg's equa- 
tion Q and ([5| change with U and V in the presence 
of the RSOC (A = 0.3) as shown in Fig. [6j where the 
data on the broken line V = -(U - 2.8) /4 in Fig. [| (U- 
V space) are exhibited. In the vicinity of the boundary 
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FIG. 6: (Color online) The eigenvalues a in the linearized 
Eliashberg's equation (|4| and Q on the broken line V = 
-(U - 2.8)/4 in Fig. |2[(C/-V space) with A = 0.3. The inset 
is a enlarged view of the region where U (V) is small (large) . 
The result is numerically obtained within the RPA. Points A 
and B correspond to those in Fig. [4] 

with the SDW or CDW phase, the eigenvalue a is large. 
Away from the boundary, the eigenvalue a is small. How- 
ever, even away from the boundary, the eigenvalue a is 
expected to increase with decreasing temperature. 

We also look at ratio between spin-singlet and - 
triplet (S z — ±1) components of the gap function on 
the broken line V = -(U - 2.8) /4 in Fig. H (U-V 
space). Fig. 7^ shows the ratio, which is defined by 
k = [A s — At / [A s + At] . Here, A s and A t denote 
the maximum absolute values of spin-singlet and -triplet 
(S z = ±1) components, respectively, k = 1 corresponds 
to purely spin-singlet pairing state while k = — 1 corre- 
sponds to purely spin-triplet one. As shown in Fig. [7] 
the ratio K slightly depends on U and V in each region, 
i.e. p + D x 2_ y 2-w&ve pairing region (0.92 < U < 1.03), 
d x 2_ y 2 + F-wave pairing one (1.03 < U < 1-48), and 
P + D xy -w&ve pairing one (1.48 < U < 2.8). Namely, 
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FIG. 7: Ratio ft between spin-singlet and -triplet (S z = ±1) 
components of the gap function on the broken line V = 
-(U - 2.8)/4 in Fig. [I] (U-V space) with A = 0.3. The 
result is numerically obtained within the RPA. Points A and 
B correspond to those in Fig. [4] 



FIG. 9: (Color online) Ratio n between spin-singlet and - 
triplet (S z = ±1) components of the gap function as a func- 
tion of A for (p+)D x 2 _ y 2-wa,ve pairing symmetry at (U, V) = 
(2.8,0), ((4,2 _ H 2+)F- wave one at (U,V) = (2.1,0.6), and 
(p+)D xv -wave one at (U, V) = (0.8,0.45). The result is nu- 
merically obtained within the RPA. 



change of pairing interaction slightly affects the ratio k 
in each region. The jump of k on the phase boundaries 
indicates that phase transitions between different pairing 
symmetries are not crossover even in the presence of the 
RSOC. 

As functions of A, the eigenvalue a in the linearized 
Eliashberg's equation ^ and {THJ) and the ratio k be- 
tween spin-singlet and -triplet (S z = ±1) components 
of the gap function are shown in Figs. [8] and [9j re- 
spectively. In both figures, data for (p+)D x 2_ y 2-w&ve 
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FIG. 8: (Color online) The eigenvalue a in the linearized 
Eliashberg's equation Q and |5| as a function of A for 
(p+)D :E 2_j ; 2-wave pairing symmetry at (U,V) = (2.8,0), 
(d a: 2_ !/ 2+)F-wave one at (U, V) = (2.1,0.6), and (p+)D xy - 
wave one at (U, V) = (0.8,0.45). The result is numerically 
obtained within the RPA. 

pairing symmetry at (U,V) = (2.8,0), (d^^+^-wave 
one at (U,V) — (2.1,0.6), and (p+)D xy -w&ve one at 
(U, V) = (0.8, 0.45) are exhibited. As shown in Fig. [8j 



the eigenvalue a slightly decreases with increasing A be- 
cause the split Fermi surface makes the nesting condition 
worse. As shown in Fig. [9j the admixture of spin-singlet 
and -triplet (S z = ±1) components is enhanced by A lin- 
early as discussed later in §111 B| 

In the above discussion, we mention the pairing sym- 
metry in the presence of the RSOC with focusing on only 
the predominant component and subcomponent whose 
amplitude is secondarily large in the gap function. Those 
always correspond to spin-singlet and -triplet (S z = ±1) 
pairings. However, we also obtain a spin-triplet (S z = 0) 
component whose amplitude is small as compared to 
those of the above components. The spin-triplet (S z = 0) 
component is always an odd function with respect to 
the Mats ubara frequency, so-called odd-frequency pair- 
ing stateP^IIHZo] Since Berezinskii proposed the odd- 
frequency pairing state in 1974,^21 it has been an impor- 
tant issue in superconductivity/superfluidity. It is ex- 
pected that the odd-frequency pairing state is discovered 
in non-centrosymmetric superconductors. 



B. Analytical discussion 

In this subsection, we discuss the pairing symmetry in 
the presence of the RSOC analytically. In the previous 
subsection, it has been clarified that pairing symmetry 
of the predominant component of the gap function in the 
presence of the RSOC is the same as in the absence of 
the RSOC. Now, we focus on pairing symmetry of the 
admixed subcomponent in the presence of the RSOC. 
We discuss below how the pairing symmetry of the sub- 
component is determined for each predominant pairing 
symmetry, i.e. spin-singlet <i 2; 2_j / 2-wave one, spin-triplet 
(S z = ±1) /-wave one, and spin-singlet d xy -~w&ve one. 

According to the linearized Eliashberg's equation (|5j), 
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if a spin-singlet component is predominant, an admixed 
spin-triplet (S z = 1) subcomponent of the anomalous 
Green's function is given by 



TT'SUb 



(k) 



Aefc [-g x (k) + ig y {k)} 



a dom 



(*)]> 

(24) 



with expanding the Green's function up to the first order 
with respect to A. Here, F™, b (fc) and A^°, m (fc) denote the 
admixed subcomponent of the anomalous Green's func- 
tion and the predominant component of the gap func- 
tion, respectively. Similarly, if a spin-triplet (S z = ±1) 
component is predominant, an admixed spin-singlet sub- 
component of the anomalous Green's function is given 
by 



TT'SUb 



(k) 



TT'SUb 



it 



(lo 2 

\ n 



r 2\2 



-g x (k)ReA^(k) 



g y (k)lmA^ m (k)] . 

(25) 



Eqs. ( 24 ) and ( 25 1 indicate that amplitude of an admixed 
subcomponent increases linearly with A. This is consis- 
tent with the numerical result within the RPA shown in 
Fig. [9} 

From the anomalous Green's function, we can derive 
the gap function by using the linearized Eliashberg's 
equation Q. For simplicity, we approximately ap- 
ply P c (k) = for c ^ 1 hereafter because the pair- 
ing symmetry is mainly determined by the term with 
P 1 (k')P 1 (—k) = 1 in the linearized Eliashberg's equa- 
tion Q in the actual calculation within the RPA. 
Then, using the convolution theorem, the linearized 
Eliashberg's equation (Wl) is rewritten as 



A^ b (fc) 



N 



F!f(r) ^F{F^(k)}, 



(26) 

(27) 
(28) 



where J 7 ' -1 ) denotes (inverse) Fourier transformation. In 
the present subsection, we focus on the pairing symmetry, 
i.e. the nodal structure of the gap function. From this 
viewpoint, it is available to neglect factors which generate 
no node in the anomalous Green's function, i.e. A/(cj„ + 



e|) 2 in Eq. Q and 4A/(w 2 1 +e|) 2 in Eq. ([25}. Thus, we 
can rewrite the anomalous Green's functions as follows, 



F% h (k) 



(—2tcosk x — 2tcos k y — fi) (sin k y 



dom 



(fc)-A ;T 



(*)], 



i sin k x ) 
(29) 



z^sub 



(*) 



TT'SUb 

r 



(— 2t cos k x 



2t cos k y — fi) 



[sinfc y ReA^ m (/c) + sinfcJmA^fc)] 



(30) 



which are effective in deriving the nodal structure of the 
gap function, by actually substituting Sk and g(k) in Eqs. 



(24 1 and |25|). Here, we neglect t' for simplicity. By using 



F^ h (k) and F^°(k) 



n 

TT'SUb 



F^ b (k) instead of F^ h (k) and 



*i (fc)— F|" (fc), respectively, in Eq. (28), we can derive 
the pairing symmetry of the admixed subcomponent of 
the gap function. 

With use of the above analytical discussion, we first 
discuss p + £>.,; 2 _ y 2 -wave pairing symmetry in the re- 
gion where U (V) is large (small). In this region, the 
dominant spin susceptibility mediates predominant spin- 
singlet c? a .2_j,2-wave pairing. We approximate the gap 
function by A^° m (fc)- A^° m (fc) = cask 



cos ky . For the 



given predominant component, we can calculate F$°{h) 
by using Eq. (|29|). For simplicity , we consider only the 

shows F^ h (r) ob- 



real part in FS^ h (k) below. Fig 



tained by Fourier transform (|28|. Note that this F|" b (r) 
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FIG. 10: 



which indicates the nodal structure of 



the admixed spin-triplet (S z = 1) subcomponent of the 



anomalous Green's function (24 1, for given predominant spin 



singlet d a ,2_ y 2-wave component approximated by A!j!° m (fc) — 
Aj° m (fc) = cosk x — cosky. Note that, strictly speaking, 
F^ b (r) displayed here corresponds to J r {ReF^ b (k)}. 



corresponds to J r {ReF|" b (A;)}. There are two kinds of 
spots in real space. One is the spot whose amplitude 
increases with \t\, which actually depends on materials 
while it is chosen as a unit of energy in pHAl The other 
is the spot whose amplitude increases with |/x|. Namely, 
the pairing symmetry of the admixed subcomponent of 
the gap function is determined by not only that of the 
predominant component but also the dispersion relation. 
On the other hand, the pairing interaction f 



the admixed spin-triplet (S z = 1) subcomponent has 
large amplitude at (x,y) = (0, 0), (±1, 0), (0, ±1) in this 
(U, V) region within the RPA. The value at (x, y) = (0, 0) 
is negative, i.e. attractive, while the values at (x,y) — 
(±1,0), (0, ±1) are positive, i.e. repulsive. Then, the re- 



sulting product f jj^(r)i 7 |" b (r) in Eq. (26) has values 



at (x,y) — (0, ±1). Sign of the value at (x, y) — (0,1) 
is opposite to that at (x, y) = (0,-1). This amplitude 
increases with \t\ at (x, y) = (0, ±1) as shown in Fig. 
Thus, the admixed spin-triplet (S z = 1) subcompo- 
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nent becomes p-wave with A|" b (/c) cx sinfcj,. Note that 
the above A|" b (fc), which is derived from the real part 
of -F™ b (fc), is real. The imaginary part of F™ h (k) gives 



the imaginary part of A™ b (fc), which is also p-wave. By 
the similar procedure, it is derived that A|" b (/c) is also p- 
wave. Thus, the admixed spin-triplet (S z = ±1) subcom- 
ponent is found to be p-wave and the amplitude increases 
with \t\. 

Next, we discuss d x 2_ y 2 + i^-wave pairing symmetry 
in the region where U and V are intermediate. In this 
region, the comparable spin and charge susceptibilities 
mediate predominant spin-triplet /-wave pairing. We 
approximate the gap function by A^° m (fc) = (cosfc x — 
cos k y ) sin k y + i(cos k x — cos k y ) sin k x . For the given pre- 



dominant component, we can calculate F*f°(k) — Ff^ h (k) 



n 
(r) 



% b (r) ob- 



by using Eq. ffl. Fig. ITT] shows F^f°(r 
tained by Fourier transform (28). On the other hand, 
the pairing interaction f jj^(r) — rjj^(r) for the ad- 
mixed spin-singlet subcomponent has large amplitude at 
(x, y) = (0, 0), (±1, 0), (0, ±1) in this (U, V) region within 
the RPA. All the values are positive, i.e. repulsive. Then, 

b (r) -Fff(r) 



the resulting product T j j.^ (r) 



77s ub 



in Eq. 

(26) has values at (x,y) = (±1,0) and (x,y) — (0,±1). 
Signs of the values at (x,y) — (±1,0) are opposite to 
those at (x,y) = (0, ±1). This amplitude increases with 
\fi\ at (x, y) = (±1, 0), (0, ±1) as shown in Fig. [TTJ Thus, 
the admixed spin-singlet subcomponent becomes d x 2_ y2 - 
wave with A*f(k) - A™ b (fc) cx cosfe* - cosk y . The 
amplitude is found to increase with \fi\. 

Finally, we discuss p + -D^-wave pairing symmetry 
in the region where U (V) is small (large). In this re- 
gion, the dominant charge susceptibility mediates pre- 
dominant spin-singlet d xy -wave pairing. We approximate 
the gap function by A^° m (fc) - Af° m (k) = sinfcj smk y . 
For the given predomin ant component, we can calculate 



F™ b (fc) by using Eq. 



(29) 



For simplicity, we consider 
Fig. 



only the real part in F|" b (fc) below 
F™ b (r) obtained by Fourier transform (28) 



this F™ b (r) corresponds to J r {Re^ b (fc)} 
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shows 
Note that 
On the other 
for the admixed 



hand, the pairing interaction F 
spin-triplet (S z — 1) subcomponent has large amplitude 
at (x,y) = (0,0),(±1,0),(0,±1) in this (U,V) region 
within the RPA. The value at (x,y) = (0,0) is negative, 
i.e. attractive, while the values at (x, y) — (±1, 0), (0, ±1) 
are positive, i.e. repulsive. Then, the resulting prod- 
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FIG. 11: F|£ b (r) - F|^ b (r), which indicates the nodal struc- 
ture of the admixed spin-singlet subcomponent of the anoma- 
lous Green's function (251, for given predominant spin-triplet 



/-wave component approximated by A!j!° m (fc) = (cosk x 
cos k y ) sin k y + i(cos k x 



(r)F|" b (r) in Eq. (26) has values at (x,y) 



(±1, 0). Sign of the value at (x^y) — (1, 0) is opposite to 
that at (x, y) = (—1,0). This amplitude increases with 
at [x, y) = (±1,0) as shown in Fig. 12 Thus, the 
admixed spin-triplet (S z = 1) subcomponent becomes 
p- wave with A|" b (fc) cx sink x . Note that the above 
A|" b (fc), which is derived from the real part of F|" b (fc), 
is real. The imaginary part of F^ h (k) gives the imagi- 
nary part of A™ b (fc), which is also p-wave. By the sim- 
ilar procedure, it is derived that A|" b (/c) is also p-wave. 
Thus, the admixed spin-triplet (S z — ±1) subcomponent 
is found to be p-wave and the amplitude increases with 

In the above analytical discussion, we derive the pair- 
ing symmetry of the admixed subcomponent of the gap 
function with use of the simplified pairing interaction, 
where some dominant modes decomposed in real space 
are chosen. The derived pairing symmetry is in agree- 
ment with that numerically calculated within the RPA in 
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state with three components 
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due to degree of freedom of spin. This relation is de- 
rived under two assumptions. One is that there is no 
inter-band pairing between the Fermi surface split by the 
RSOC. The other is that intra-band pairings on the in- 
ner and outer Fermi surface have same pairing symme- 
try. In the present calculation within the RPA, the for- 
mer assumption is broken, i.e. there is finite inter-band 
pairing between the Fermi surface split by the RSOC. 
Therefore, the previous discussion (31 1 is not applica- 
ble. Actually, the previous and present studies give dif- 
ferent conclusions. According to the previous study, spin- 
triplet /-wave pairing symmetry should be admixed with 
spin-singlet d x 2_ y 2-wecve one. On the other hand, in the 
present study on the basis of the RPA, it is derived that 
spin-triplet (S z = ±1) p-wave pairing symmetry is ad- 
mixed with spin-singlet d x 2_ y 2-wave one in the region 
where U (V) is large (small) as shown in Fig. [4] Thus, 
in the case where the above two assumptions do not hold, 
we must take into consideration the dispersion relation 
and the pairing interaction. 



FIG. 12: F-ff? (r), which indicates the nodal structure of the 
admixed spin-triplet (S z = 1) subcomponent of the anoma- 
lous Green's function (241, for given predominant spin-singlet 



d^-wave component approximated by A5? m (fc) — A?2 m (fc) 



sin k x sin k y . Note that, strictly speaking, F^ b (r) displayed 
here corresponds to J r {ReF|^ b (fc)}. 



IV. SUMMARY 

In order to study the pairing symmetry in non- 
centrosymmetric superconductors, we have solved the 
linearized Fliashberg's equation on the two-dimensional 
extended Hubbard model in the presence of the RSOC 
within the RPA. We found that three types of the pair- 
ing symmetries appeared in the U-V phase diagram in 
the presence of the RSOC. In the region where U (V) is 
large (small), p + D x 2 _ y 2 -wave pairing symmetry, which 
is predominantly spin-singlet (4,2 _ y 2 -wave one admixed 
with spin-triplet (S z = ±1) p-wave one, is the most sta- 
ble. In the region where U and V are intermediate, 
d x 2_ y 2 + F-wecve pairing symmetry, which is predomi- 
nantly spin-triplet (S z = ±1) /-wave one admixed with 
spin-singlet d x 2_ y 2-wave one, is the most stable. In the 
region where U (V) is small (large), p + D xy -wsve pairing 
symmetry, which is predominantly spin-singlet d xy -wsve 
one admixed with spin-triplet (S z = ±1) p-wave one, is 
the most stable. 

From analytical study, we found that pairing symme- 
try of an admixed subcomponent of the gap function de- 
pends on not only that of the predominant component 
but also the dispersion relation and momentum/space 
dependence of the pairing interaction. Amplitude of the 
admixed subcomponent of the gap function depends on 
the dispersion relation, i.e. the hopping t and the chem- 
where d{k) is the gap function for spin-triplet pairing ical potential [i, as follows. For the p + D x 2_ y 2-wave 



terms of the symmetrical class such as p- or d-wave sym- 
metry. However, there is difference in the nodal structure 
in detail. The difference comes from the simplification of 
the pairing interaction. Considering the higher order har- 
monic components in the pairing interaction, of course, 
we can reproduce in detail the nodal structure in the ad- 
mixed subcomponent of the gap function obtained by the 
numerical calculation within the RPA. 

There exists previous analytical discussion on the ad- 
mixture of the pairing symmetryP According to the pre- 
vious study, spin-triplet components are related to spin- 
singlet one in the gap function by 



#)«^|[A ti (fc)-A ;t (fc)], (31) 
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pairing symmetry, amplitude of the admixed spin-triplet 
(S z = ±1) p-wave subcomponent increases with \t\. For 
the d x 2_ y 2 + F-wave pairing symmetry, amplitude of 
the admixed spin-singlet d x 2_ y 2-wave subcomponent in- 
creases with For the p + D xy -w&ve pairing symmetry, 
amplitude of the admixed spin-triplet (S z = ±1) p-wave 
subcomponent increases with \fj,\. 
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